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Abstract 

Diagonalizable pseudo-Hermitian Hamiltonians with real and discrete spec- 
tra, which are superpartners of Hermitian Hamiltonians, must be r^-pseudo- 
Hermitian with Hermitian and invertible rj operators. We show that de- 
spite the fact that an rj operator produced by a supersymmetric transfor- 
mation, corresponding to the exact supersymmetry, is not invertible in the 
whole Hilbert space, it can be used to find the eigenf unctions of a Hermi- 
tian operator equivalent to the given pseudo-Hermitian Hamiltonian. Once 
the eigenfunctions of the Hermitian operator are found the operator may be 
reconstructed with the help of the spectral decomposition. 
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After works by Andrianov and Cannata with coworkers the supersym- 
metric quantum mechanics (SUSY QM) became an effective tool in study- 
ing different properties of the Schrodinger equation with a complex-valued 
potential. Thus in j2[ the author shows that new exactly solvable complex- 
valued potentials may be obtained with the help of second order irreducible 
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SUSY transformations. In a possibility is discussed to transform non- 
diagonalizable Hamiltonians into diagonalizable ones and vice versa and in 
the author stresses that SUSY transformations may create Hamiltonians 
with spectral singularities inside the continuous spectrum. Green functions 
for some complex SUSY partners of real potentials are calculated in [sj . In 
SUSY QM is applied to study a dynamo effect and in j3] spontaneous break- 
down of PT symmetry is related with the presence of the supersymmetry in 
the system. 

One of the main features of supersymmetric transformations correspond- 
ing to the exact supersymmetry is that the transformation operator has a 
nontrivial kernel. This results in the property that a non-Hermitian di- 
agonalizable Hamiltonian with a real and discrete spectrum, which is a 
super-partner of an initial Hermitian Hamiltonian, becomes (weak) r^-pseudo- 
Hermitian where 77 = r/Ms not invertible in the whole Hilbert space. The 
possibility for i] to be non-invertible was first considered by Fitio |8| and lat- 
ter explored in [oj in the context of generating non-Hermitian Hamiltonians 
with real spectra. Nevertheless, these authors completely skipped the ques- 
tion whether or not a new Hilbert space, where the non-Hermitian operator 
becomes Hermitian, may be constructed with the help of such an rj operator. 
Indeed, if 77 = 0^0 has a nontrivial kernel, kero 7^ 0, a new inner product 



defined as (see e.g. |10| ) 



((V^,0)):=(V^|r7|0), v = V^ (1) 

becomes positive-semidefinite even if the Hamiltonian is diagonalizable. In 
view of a theorem proved in ll| this means that such an t] operator is 
not suitable for constructing the new inner product and there should exist 
another rj operator which is invertible in the whole Hilbert space. 

The main aim of the current Letter is to show that such an rj operator 
may be used for constructing a Hermitian operator equivalent to the given 
non-Hermitian one acting in the same Hilbert space. Thus, the approach 
based on equivalence classes (see lill proves to be more advantageous in 



this case than the usual scheme [10|, lUj based on Eq. ((T 



For simplicity we will consider a regular Sturm-Liouville problem defined 
by a differential expression h = —d1 + Vo(x) {dx = ^) with a real-valued 
potential function Vo(x), x G [0,d] and the Dirichlet boundary conditions at 
X = 0,d. Under rather mild conditions on the function Vo{x) {Vq{x) is real 
and has a bounded derivative, see e.g. [isj) this differential equation together 
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with the Dirichlet boundary conditions defines a selfadjoint operator in the 
Hilbert space L'^{0,d) of square integrable functions over the interval [0,d]. 
To simphfy notations we will denote this operator by the same symbol h = . 
Then the corresponding boundary value problem 

h%l) = E%l), ^/>(0) = = (2) 

(and, hence, operator h) has an unbounded real and discrete spectrum [l^ 
E = En = E* with eigenfunctions ip = ipn = i'n, = 1, 2, . . . 

Let Dh be the domain of definition of h. It consists of the twice con- 
tinuously differentiable functions defined Vx G (0, d) satisfying the Dirichlet 
boundary conditions. Let L be a first order differential intertwining operator 
between h and a non-selfadjoint operator H = —d^ + V{x) 

Lh = HL (3) 

where in general V{x) is a complex-valued function which we will describe 
below. Actually, for the moment we consider H only as a differential expres- 
sion but below we will supply it with proper boundary conditions which will 
transform H into a non-selfadjoint operator acting in the same Hilbert space 
L'^{0,d). Operator L is defined in terms of a complex- valued superpotential 
W = W{x) as follows L = —d^ + W where W = u^/u {ux := du{x)/dx) with 
u{x) being an (essentially) complex (i.e. W must be complex) and nodeless 
solution to the differential equation hu = au, a G M and Lu = (see e.g. [l|). 
Note that if a coincides with one of the levels E^, n = 1,2, . . ., i.e. a = E^q, 
the Hamiltonian H becomes non-diagonalizable (see e.g. |3]), a possibility 
we would like to avoid. Therefore we will assume a ^ E^, n = 1,2,.... 
If now ip = ip{x) is a solution to the differential equation {h — E)ip = 
then (f = ip{x) = Lip{x) satisfies differential equation {H — E)(p = 0, 
H = —dl + V{x), V{x) = Vo(x) — 2Wx- By this reason we have 

ip, + W^={E-Vo + W, + W^)ij . (4) 

From (jl]) it follows that if ip satisfies the Dirichlet boundary conditions (see 
([2D) then 

{ipx + W^) (0) = + W^) (rf) = . (5) 

From here we conclude that the functions = Lipni n = 1,2,... solve the 
following boundary value problem 

H(fn = Enifn , {Vnx + W ipn){Q) = {^nx + WLPn){d) = . (6) 
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This means that in the Hilbert space L^{0,d) a non-selfadjoint operator H 
with the domain of definition Dfj is defined. Domain consists of the twice 
continuously differentiable functions defined Vx G (0, d) satisfying boundary 
conditions ([5]). It is important to note that the function ipo = 1/u also solves 
the boundary value problem ([6]). This means that the spectrum of H consists 
of the levels n = 1, 2, . . . plus E = Eq = a ^ E^- From the point of view 
of SUSY QM this corresponds to the exact supersymmetry (see e.g. [U]). 

In what follows we do not need to use the full domains of definition of 
the operators h and H. Instead we will use domains T>h and T>h which are 
the sets of all finite linear combinations of the eigenfunctions of h and H 
respectively. Operator L is well defined on and maps — )• Vh. 

Operator W Hermitian adjoint to H is defined by the adjoint boundary 
value problem 

H^^ = Ei , (e. + W*0 (0) = (e. + W*0 (d) = (7) 

with W = —d1 + V*{x). Its domain of definition D^^ consists of the twice 
continuously differentiable functions defined Vx G (0, d) with the boundary 
conditions as given in ([7]). The spectrum of H'^ coincides with the spectrum 
of H, E = En, n = 0,1 . . .. Its eigenfunctions ^ = C,n = V^n, n = 0, 1, . . . form 
a basis in L^{0,d) biorthogonal to {(pn} 

(^nlv^m) = 0, n 7^ m. (8) 

In what follows we will use Vjjt which is the set of all finite linear combina- 
tions of iJ^ eigenfunctions ^n- 
From (|3]) it follows that 

hL^ = L^H^ (9) 

where L'^ = 8^ + W*. This relation means that transforms eigenfunctions 
of if^ to the eigenfunctions of h except for = V^o ~ l/'u* which enters the 
kernel of L\ L'^C.o = 0. Operator L''" is well defined on P^/t. It realizes the 
mapping V^t Vh- 

From (E]) and (jH]) one deduces that 

HLL^ = WH^ . (10) 

Since the spectrum of H is real, it should be r^-pseudo-Hermitian llj, i.e. 



there should exist an invertible Hermitian operator rj such that if ^ = rjHrj ^ 
or equivalently 

r]H = H^r]. (11) 
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If LL'^ were invertible, comparing (ITU]) and fllip we could conclude that r] = 
{LL^)~^. Unfortunately this is not so. Nevertheless, as we show below, the 
operator LL^ can be used to calculate eigenfunctions of a Hermitian operator 
ho = hi equivalent to H'^ (and, hence, to H). Then ho may be reconstructed 
with the help of the spectral decomposition. Note that composition LL^ is 
well defined on Vj^t and maps Vjjt ^ Vh while its domain of definition is 

Operator LL^ is a second order differential operator, 

LL^ = -dl + {W - W*)d^ + WW* - W* (12) 

and in the space V^j] it has a one- dimensional kernel with the basis function 
= 1/u*. From f IlOp it follows that it transforms eigenfunctions of H'^ 
to eigenfunctions of H except for E = a. Therefore using biorthogonality 
relation ([8]) one gets 

(a|LLtU = 0, n^m. (13) 

It is not difficult to see that being defined on D^^ operator LL^ is selfadjoint 
and positive semidefinite. Therefore the boundary value problem 

LLtH = A^H , (S, + W-E) (0) = (S, + W*'E) (d) =0 (14) 

has a real non-negative spectrum = A^, n = 0, 1, . . ., Aq = with eigen- 
functions S = En, n = 0, 1, . . ., Ho = (Po- Its Hermitian square root is well 
defined in L^{0, d) and may be written using its spectral decomposition 

oo 

(LLt)i/2 = ^A,|Hfc)(Sfc|. (15) 

Let us consider the functions 

$„ = (LLt)i/2^„, n = l,2,..., $o = So. (16) 

From (fT3|) it follows that =0, n ^ m, n,m = 1,2. . .. Moreover, 

($„|$o) = 0, n = 1, 2 . . . since LL^Eo = 0. Thus the functions form 
an orthogonal set. Furthermore, since n = 0, 1, . . . form a basis in -L^(0, d) 
the functions f lT6|) form an orthogonal basis in L^(0, d). 

Denote T>hoi the space of all finite linear combinations of the functions 
n = 1, 2, . . .. Being restricted to this space, {LL'^Y^'^ becomes invertible. 
Therefore operator 

hoi = {Wf/^H^LL^)-^/^ (17) 
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is well defined on 2^/ioi- Moreover, from intertwining relation fllOp it follows 
that it is Hermitian on "Ph,,!, /iqi = ^oi and /loi'^'n = En^n, n = 1,2,.... 
Therefore Eq. f|T7|) may be rewritten as a couple of intertwining relations 

(LL^y^'^hoi = H{Wf'^ , {Wf'^H^ = hoiiLL^y/' . (18) 

Eqs. f|T8l) mean that the operator (LL"!")^/^, being restricted to Vh^j^, real- 
izes an equivalence transformation between a restriction of H and /iqi- More 
precisely it transforms eigenfunctions $„ of /iqi to eigenfunctions ipn of if, 
{LL'^Y^'^^n = '■Pn, n = 1,2,.... Similarly, since (LL'^Y^'^ transforms eigen- 
functions of H^^ to eigenfunctions $„ of /iqi, (LL'^Y^'^^n = '^'n, n = 1, 2, . . ., 
it realizes an equivalence transformation between hoi and a restriction of if^. 

Now we will extend operator hoi from the space Vho^ to the space P/ig of all 
finite linear combinations of the functions ^n, n = 0, 1, . . .. By construction 
one has = T^hoi © T^Hqo where P/^jq is a one-dimensional space with the 
basis function $o which is orthogonal to T^h^^. By this reason if we fix the 
additional eigenvalue of ho to be equal to -Eq = a, this continuation, if 
restricted to be selfadjoint, becomes unique 

oo 

/io = 5^|iV„|2^„|$„)($„|. (19) 

n=0 

Here Nn is a normalization coefficient of the function $„. Although this 
approach permitted us to construct a unique operator ho equivalent to H\ 
an invertible t] operator may lead to another h^ which is unitary equivalent 
to (dH). 

Let us consider the simplest case when the boundary value problem 
(fT4|) may be solved exactly. It corresponds to transformation function u = 
exp (iax) with W = ia, a = a^. The initial boundary value problem ([2]) 
in this case has zero potential. For the transformation operator one gets 
L = —dx + ia and LL^ = —d^ + 2iadx-\- a? ■ The spectrum and eigenfunctions 
of LL^ are 

Xl = kl, kn = iTn/d, n = 0,1,2,... (20) 

cos{knx), n=l,2,..., Ho = y^e'"'^. 

Since W does not depend on x the operator H has still zero potential 
but it is non- selfadjoint because of boundary conditions ([5]) which in this 
particular case read 

{ipx + iav) (0) = + ia^) (d) = . (21) 
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We recognize here the boundary value problem generated by the Robin 
boundary conditions previously studied in [l5|, ll6| . The spectrum of H con- 
tains all energy levels corresponding to the boundary value problem (|2]) with 
zero potential En = k^, n = 1,2,... plus the level Eq = a = and the 
supersymmetry is exact. The corresponding eigenfunctions (not normalized 
to unity) read 



ta 



ifn = cos{knx) - — sin(A;„x) , n = 1, 2, . . . , ipo 



Eigenfunctions of if^ are simply complex conjugated eigenfunctions of H, 
= ^*n- As it was mentioned above if a = E^, n = 1, 2, . . . both Hamiltonian 
H and H'^ become non-diagonalizable and therefore we have to assume a ^ 
TTn/d, n = 1,2, . . .. 

To find eigenfunctions $„, n = 1,2,... according to (ITB]) and f ll6p we 
have to evaluate integrals Smn = (Sml^n), which are almost standard (see 
e.g. [13) 



2iadV2de~''"^{{-l)'^+'' - e''"^)(a^d^ - mV^) 



(22) 



and then to sum up the series (fT5|) with (LL'^Y^'^ applied to Using Eqs. 
(12 op and (I22p one can reduce this series to a combination of two known series 
(see ^12]) 

— ; — = /3(p) cos[(7r — x)p] + - / cos[(p — l/2)t — xp] CSC -(it, (23) 



n=0 



n + p 



n=0 



^— 1)" COS nx 
n + p 



/3(p) cos[xp] - - / sin[(p - l/2)t - xp] sec -(it . (24) 
2 Jo 2 



Finally for the eigenfunctions of we get 



Aad 



TT 



2nd7i 
7(C) cos 



XTC( 



— 7((5) cos 



7r5 



X7[6 



7 



+ 



+ 



+ 



TT 
TT 
TT 



COS 



7(5) COS 



TTX6 



cos 



Here n = 1, 2 . . ., 5 



ad 



COS 



7(0 COS 



TlxS 
7TX( 



t6 



cot -(it 



t(5 



cos 



d 



tan -dt 
2 



ad 



+ n and 



1 /rf^ 



2 

2+11 



m 



with r[z] being the Euler gamma function. 

To summarize, we have shown that despite the fact that the operator 
LL^ produced by SUSY QM with exact supersymmetry is not invertible, 
it can be used to find eigenfunctions of Hermitian operator equivalent to 
non-Hermitian operator generated by SUSY QM. Once these eigenfunctions 
are found the Hermitian operator may be reconstructed using its spectral 
decomposition. 
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